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Abstract 

We explore the pseudoscalar rji, and the scalar Xbo decays into to probe whether it is possible 
to probe the Higgs sectors beyond that of the Standard Model. We, in particular, focus on the 
Minimal Supersymmetric Standard Model, and determine the effects of its Higgs bosons on the 
aforementioned bottomonium decays into lepton pairs. We find that the dileptonic branchings of 
the bottomonia can be sizeable for a relatively light Higgs sector. 
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I. INTRODUCTION 



Having the LHC started, the search for physics at the terascale has entered a new phase. 
The ATLAS and CMS experiments at the LHC will search for new particles and forces 
while LHCb will provide a more accurate description of flavor physics. Each experiment, 
combined with others, will provide important information about nature of new physics 
awaiting discovery. It is thus rather timely to discuss and analyze ways of extracting TeV 
scale physics in relation to the measurements at the LHC experiments. 

In terms of its content and goal, the present work falls in the interface between flavor 
physics and Higgs physics in that we aim at exploring finger prints of yet-to-be discovered 
Higgs sector (to be discovered at the CMS and ATLAS experiments in the leptonic decay 
distributions of heavy hadrons (to be accurately measured at the LHCb experiment j2|). 

At present, we do not have any clue of what Higgs sector is awaiting for discovery at 
the LHC. On the other hand, the experiments at B factories have, by now, established a 
grand view of the flavor physics. The experimental precision is increasing steadily and has 
already started challenging our understanding of the flavor violation. Over the years, various 
B meson decay rates and charge asymmetries have been measured and novel quarkonium 
states have been discovered. The B meson inventory of the existing storage rings comes from 
the decays of {bb) states (bottomonium states) produced at asymmetric electron-positron 
collisions {e.g. PEPII at SLAC and KEK-B at KEK). Of course, all kinds of bottomonia 
with varying spin and CP quantum numbers will be produced at the LHC in gluon-gluon 
or gluon-gluon-gluon fusion channels. 

In principle, one ought to use every single opportunity to extract information about other 
sectors of a given theory by using the available information from B physics. Examples of 
such efforts involve quark EDMs and flavor- violation Higgs connection [4]. The radia- 
tive, leptonic or semileptonic decays of hadrons are particularly suitable for strengthening 
experimental identification and theoretical prediction, and thus, in this work we attempt 
at answering the following question: By measuring the decay rates of certain (bb) 
states, preferably but not necessarily into , can we establish the existence 
and nature of Higgs bosons? The choice of bottomonium system stems from not only 
its perturbative nature but also its appreciable coupling to Higgs fields. 

In what follows, in regard to the question raised above, we will study a generic Higgs 
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sector extending that of the SM. In the next section, we will provide an explicit discussion of 
the dileptonic Bottomonium decays into lepton pairs. In Sec. 3 we will numerically analyze 
the decay rates by taking MSSM to be the new physics candidate model and fixing the 
unknown parameters to two different data sets taken from LEP indications and from SPSla 
point. In Sec. 4 we conclude. 



II. DILEPTONIC BOTTOMONIUM DECAYS 



The quarkonium systems have been under intense study since the discovery of the charm 
quark ^. That light MeV-mass Higgs bosons could be produced in quarkonium decays was 
first discussed in [g, [?]. The decays of additional TeV-mass heavy quark bound states into 
fermions as well as Higgs and gauge bosons have been analyzed in [Sj. In this work we will 
discuss Higgs boson search via dileptonic bottomonium decays. The two sides, hadronic and 
Higgs aspects, of our discussions can be described as follows: 

1. We will focus on bottomonium states, in particular, the pseudoscalar rji) (an 5*- wave 
J^'^ = 0~"^ state) and the scalar Xbo P-wave J^'" = 0~^~^ state). Unlike the char- 
monium system where such states have already been experimentally established, the 



experimental efforts still continue to estab. 
they have already been observed 0.0' 



ish quantum numbers of rjh and Xbo though 
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121]. The experiments at Tevatron, B 



factories and LHC are expected to fully construct and measure partial widths of these 
b b mesons. 

2. We will not restrict ourselves to standard model Higgs sector. In fact, we will consider 
models with two Higgs doublets, and Hfi one giving mass to down-type fermions 
other to up-type fermions, as encountered in the MSSM. (One may, of course, consider 
more general Yukawa structures |lj].) The spectrum consists of three Higgs bosons: 
the CP-even ones h and H and a CP-odd one A. Their interactions with b quark and 
charged leptons are given by 

= gl7fh + g{jJfH + g{Jt^,fA (1) 
where f = b,i, and the Yukawa couplings Qx are given by 

gl = hf'^^ [sin(/3 - a) - tan/5 cos(/5 - a)] 



gf^ = /i^^ [cos(/9 -a)+ tan /3 sin(/5 - a)] 

g{ = hf' tan p. (2) 

Here tan/? = {H^)/{H^), a is the mixing between - and - (Z/^) such 
that a = 1/2 aicsm[-{m\ + ml) /{mjj- ml) sm2P], h^'^^ = {g2mf)/{2Mw) is the 
Yukawa couphng of fermion / in the SM. If there exists exphcit CP violation sources 
in the theory then none of the Higgs bosons can possess definite CP quantum number, 
and thus they couple to fermions as f{a + ilry^)f as in, for instance, the MSSM with 
complex soft terms with one-loop Higgs potential 

Having specified the framework in both Higgs and meson sides, we now turn to an explicit 
computation of the decay rates of bottomonia. In this respect, the decay rates of r^b and xw 
into lepton pairs are then given by 




+ 4v/^— ^^^^^^ 



(l-rA)(l-rz) 



M4 

X60 



(3) 



\l-rh I-thJ 

where -Rs(O) (i?p(0)) is the S-wave (derivative of P-wave) quarkonium wavefunction at 
the origin [Sj, gz = e/(4sin6'vy cos6'vy), rj = mj/Mjr {i = i, A, h, H, Z), and Pi = 
(1 - 4m2/M2 )^/^ where X = r]^ for r]^ -> £+£- and X = Xbo for Xbo ^ 
From these decay rates one notes that: 

1. Thanks to their J^'-' structures, the two bottomonia, rjh and Xbo, explicitly distinguish 
between CP=+1 and CP=-1 Higgs bosons. This aspect proves very important for 
establishing the nature of the Higgs bosons as well as structure of the non-SM Higgs 
sector at the LHC and its successor NLC (see [sl for a detailed discussion of different 

jPC 

mesons . 

2. The couplings of the Higgs bosons to down- type fermions experience big enhancements 
at large tan (3 as preferred by LEP experiments. Indeed, contributions of A and H 
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grow as (tan/3/Mj7^^)^ which can provide a detectable signal for colhder experiments 
such as the LHCb. 



3. As is seen from ([3]), as a direct consequence of the quantum numbers of rjh meson, the 
decay rjf, exclusively involves the vector bosons and pseudoscalar Higgs bosons. 
On the other hand, again due to its quantum numbers, the decay Xbo ~^ singles 
out the CP=+1 Higgs bosons. (Nevertheless, one keeps in mind that Xbo ~^ can 
exhibit a non- negligible dependence on the masses of the pseudoscalar Higgs bosons 
depending on the details of the mixing angle a of the CP-even Higgs sector). This 
prime difference between the two mesons proves highly useful for probing the nature of 
the 'new physics' that will be discovered at the LHC. Depending on the nature of the 
deviations from the SM expectations, one might determine, within the experimental 
uncertainities, whether the new physics involve new pseudoscalar Higgs bosons (like 
the MSSM or NMSSM) or new scalar Higgs bosons (like MSSM or U(l)' models or 
NMSSM) or new gauge bosons (like U(l)' invariance of left-right symmetric models). 

4. In (l3l) we have focused particularly on extended Higgs sectors (taken to be a generic 
two-doublet model fitting to the Higgs sector of the MSSM). However, one can consider 
extended gauge sectors as well. In this case, similar to the ZjW boson contributions, 
one expects anomalous behavior in — * i^l" (compared to the SM prediction) to 
arise also from extended gauge sectors containing Z' jW' gauge bosons. In this work 
we will not investigate this option sice experimental bounds force Z' jW to stay heavy 
(though in realistic models the Higgs sector itself behaves differently isj]). 

5. In the decoupling limit , it turns out that /? — a~7r/2in which case h behaves 
as in the SM yet H and A Higgs bosons possess tan /3-enhanced Yukawa interactions. 

In the next section we will perform a numerical study of the decay rates ([3]) in view of 
disentangling H and A effects from the rest. The analysis, once confirmed experimentally, 
might provide important information about the nature of the Higgs sector awaiting discovery 
at the LHC. 
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III. NUMERICAL ANALYSIS 



In this section we analyze the decay widths discussed above numerically. In doing this, 
the SM prediction for the decay rate will be compared with those of the MSSM for Xbo and 
rjj, decays, comparatively. In particular, we take Higgs boson of the SM degenerate in mass 
with the lightest Higgs boson of the MSSM, and consider the ratios 

and 

(5) 

in making the numerical estimates. The parameter values for which these ratios exceed 
unity significantly are expected to yield observable signals. In course of the analysis we 
scan the parameter space of the MSSM Higgs sector by varying rrih, rriH, mA and tan/3 in 
a considerably wide range. We focus on two parameter ranges: 

• SUSY Parameter Space I (SPSI): 

mh = 98 ±5 GeV, mn = 115 ± 5 GeV, mA = 89 ± 5 GeV, tan p = 10 ± 2.5 (6) 



which is inspired from the reanalysis of the LEP results mentioned in (l5i |. 
• SUSY Parameter Space II (SPSII): 

m;, = 1 15 ± 5 GeV, mn = 425 ± 5 GeV, mA = 424.9 ± 5 GeV, tan /? = 10 ± 2.5 (7) 
which is inspired from the SPSla parameter space of the MSSM. 

In plotting a particular figure we vary one parameter while keeping the rest at their mid- 
values. Depicted in Fig. 1 (for SPSI) and Fig. 2 (for SPSII) are the ratios in (j4]) and ([5]) as 
a function of the lightest Higgs boson mass m/^. As can be seen from the left panels of the 
figures, the ratio of the SUSY prediction to the SM prediction does not vary for the rjb decay 
and these ratios are approximately 1.23 and 1.02 for the Figs. [T]and[2l respectively. This can 
be easily understood from (jl]), to which the CP-even Higgs bosons do not contribute at all. 
The Tjb decay would probe new CP=-1 Higgs bosons and new gauge bosons, as can be seen 
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FIG. 1: Variation of the decay rate ratios against the Ughtest CP-even Higgs mass mh for r]b (left 
panel) and Xbo (right panel), for the SPSI parameter space. 

from the same equation. Nevertheless, the difference persistent in the MSSM's prediction 
can be an important clue for the future measurements. 

While the impact of different data sets are sensible for the rji) decay, for the xto decay it 
turns out to be much stronger as can be seen from the right panels of the same figures. For 
instance the impact of increasing the mass of the lightest CP-even Higgs boson can enhance 
the SUSY/SM ratio from 5050 to 7700 for SPSI parameter set. Similarly, it increases from 
62.5 to 84.5 for SPSII set. Here, as a result of the quantum numbers of Xbo, the contributions 
of the h and H bosons become clearly visible, which can be seen from ([5]). Normally, as 
mh increases the Vh, and hence, the decay rate ratios increase but the dominant behavior is 
determined by the coupling terms. In any case, the prediction of the MSSM is very large 
than that of the SM prediction, which makes the Xbo decay a promising candidate for probing 
the new CP-even Higgs bosons of the 'new physics'. 

Depicted in Figs. [3] and HI are variations of the decay rate ratios with the heavy CP-even 
Higgs boson mass, tuh- The general behavior is similar to those in Figs. [1] and [21 except 
that the Xbo decay ratio decreases as mn increases, for both of the parameter sets. This can 
be understand from ([5]) wherein the decay rate ratio is inversely proportional to m//^. 

Depicted in Figs. [5]and[6]are the ratios of the decay rates as a function of the pseudoscalar 
Higgs boson mass rriA- As suggested by the left panels of the figures be (the rjb decays) the 
ratios decrease with increasing niA to a small extend, for both of the parameter sets. For the 
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FIG. 2: The same as FiglU but for SPSII parameter set. 
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FIG. 3: Variation of the decay rate ratios against the heavy CP-even Higgs mass mn for % (left 
panel) and Xbo (right panel), for the SPSI parameter space. 

Xbo decays, the reaction response to variation in rriA is much more pronounced: The decay 
rate ratio ranges from 6325 to 6760 for SPSI, and does from 71.85 to 73.70 for SPSII. It is 
important to stress that, as suggested by formulae ([5]), the dependence of the Xbo decay 
follows from the dependencies of the H and h couplings on the Higgs mixing angle a. 

Another important parameter for the MSSM is the ratio of the vacuum expectation values 
of the Higgs bosons, the tan/5. The tan/5 dependencies of the related decays are depicted 
in Figs. [7] and M for rji, and Xbo for the two parameter sets employed in previous figures. 
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FIG. 4: The same as FigjSl but for SPSII parameter set. 
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FIG. 5: Variation of the decay rate ratios against the CP-odd Higgs boson mass for % (left 
panel) and Xfeo (right panel), for the SPSI parameter space. 

For the tan/3 values contained in SPSI and SPSII, r/b decay exhibits less sensitivity to 
tan [3 compared to xw decay. Nevertheless, the r]i, decay stands as a sensitive probe of tan /? 
since it scales as (tan /5)^, which becomes quite sizeable at large values of tan j3. For instance, 
as can be seen from the left panel of Fig. [7| MSSM's prediction can be ~ 1.5 times larger 
than of the SM for reasonable values of tan (3. The impact of the tan (3 variable for the Xbo 
decays is always supportive to claim that the MSSM prediction can be four (right panel of 
Fig. [7]) or two orders (right panel of Fig. [H]) of magnitude larger than the SM results. 
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FIG. 6: The same as FigjSl but for SPSII parameter set. 
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FIG. 7: Variation of the decay rate ratios against tan/3 for % (left panel) and Xbo (right panel), 
for the SPSI parameter space. 

Our last figure is devoted to examining the decay rates in the decoupling limit i.e. the 
domain in which itia = fnn and it is much larger than rrih. For this aim, we take tan f3 = 10, 
rrih = 120 GeV and vary rriH = itla from 124 to 920 GeV. The numerical results are 
depicted in Fig. [9l As can be seen from the left panel of the very figure, rjf, ratio decreases 
as tuh = friA increases and its prediction does not offer a difference more than ~ 3%. The 
largest effect occurs when mn = rriA is not much larger than rrih, which actually means that 
the rjf, decay cannot give any significant result in the decoupling regime. 
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FIG. 9: Variations of the decay rate ratios against uih = tha. Here we fix the parameters as 
tan /? = 10 and rrih = 120 GeV. 

Coming to Xbo, however, one notes from the right panel of Fig. [9] that, Xbo is very 
sensitive to the variation of rriA = fnn, especially for low values of the heavy Higgs bosons. 
As rriH converges to rrih the ratio F^^^^ [xbo ~^ i^i^) /F^^ {xbo — ^ can be enhanced 

up to ~ 900. Of course this ratio can be further enhanced by increasing the tan/3. The 
lesson from this figure is that the Xbo decay in supersymmetry is a candidate with 

significantly enhanced predictions with respect to the standard model rate. 
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SM 




MSSMfSPSI) 


MSSMfSPSIIl 


% ^ £+t 


Cornell 


5.08 X 10~ 


-7 


6.23 X IQ-'^ 


5.13 X lO^"^ 




Richardson 


2.37 X 10" 


-7 


2.91 X IQ-'^ 


2.39 X IQ-'^ 




Wisconsin 


1.86 X 10" 


-7 


2.29 X 10"^ 


1.88 X IQ-^ 




Coulomb 


1.02 X 10" 


-7 


1.25 X lO"'^ 


1.03 X 10-7 



TABLE I: The branching ratios of the rjb decay for different potential models [8i] in the SM and in 
the MSSM. 



Decay 


Potential 


SM(SPSI, SPSII) 


MSSM(SPSI, SPSII) 


Xbo ^ i^i- 


Coulomb 


(6.20 X 10"^^, 3.25 X 10"^^) 


(4.03 X 10-^2^ 2.36 X 10"^^) 




others 


(6.20 X lO-^'', 3.25 X 10"^^) 


(4.03 X 10"^°, 2.36 X 10"^^) 



TABLE II: The branching ratios of the Xbo decay as in Tab. HI 

Moreover, we estimated the branching ratios of the rjb and Xbo decays into i^i" pairs for 
both SM and MSSM processes using two parameter spaces: SPSI and SPSIL In doing this, 
different potential model wave functions are examined Q| to probe the arbitrariness in the 
potential dependency. Our findings are presented in Tabs. [T] and [III for m and Xbo decays. 



respectively. As input parameters we used F^^^ = 10 MeV taken from [l6| and T^^^ = 320 



keV from 



17|. 



As can be read from Tab. [T] our predictions for the branching ratios of the rjb 
decay in the MSSM is ~ 1 (SPSII) and 1.2 (SPSI) times larger than the SM values. On the 
other hand, as can be read from Tab. [TTl MSSM predictions for the Xbo ~^ branching 
ratio can be as large as 73 (SPSII) or even 6500 (SPSI) times larger than that of the SM 
predictions . 

It should be noticed for both of the decays that they are rare decays. It is possible to 
enhance the related predictions theoretically in the MSSM, as examined in this section, but 
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experimental verification of sucli predictions is a cliallenging task. 



IV. CONCLUSION 

In tliis work we liave studied dileptonic bottomonium decays in regard to tlieir sensitivity 
to Higgs bosons of either CP quantum number. We have found that, dileptonic branching 
of Xbo is a highly sensitive probe of the extended Higgs sector in that the rate increases 
significantly compared to the SM prediction. 

Theoretically, comparison of the rih ratio with the Xbo ratio shows that, for the selected 
parameter ranges, the likelihood of observing the Higgs bosons via dileptonic rji, decays turns 
out to be much smaller than those of Xbo decays. On the other, experimentally, since the 
predictions of the branching ratios are at the order of ~ 10^^ for the rj and ~ 10~^° for the 
X decays, both in the SM and in the MSSM, r^f, turns out to be a better candidate for the 
observation of the Higgs bosons over these rare decays. 

The results found here, given the high-luminosity, high-energy nature of the LHC exper- 
iments, can be tested at the LHCb experiments is not at the B factories. Such a test, if 
conducted, would provide a confirmation strategy if not a discovery strategy for extended 



Higgs sectors. The recent paper by [l8| also discusses the bottomonium decays with partic- 



ular emphasis on light pseudoscalar Higgs which can be realized in the NMSSM. 
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